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1. Introduction 

It can be argued that the strong coupling constant a s , or, equivalently, the QCD scale Aqcd 
are the fundamental quantities which require a high level of precision as they enter in most cal- 
culations, both perturbative and non-perturbative. They can be extracted both from experiment or 
numerically, with the help of Lattice QCD. The latter is a complicated procedure and it is imperative 
to understand the underlying mechanics to estimate correctly the systematic errors. For the review 
of such calculations see, for instance, [1, 2, 3, 4, 5, 6, 7]. We presented our most recent results 
in [8] and here we will demonstrate that our approach and estimation of systematics is justified. 

The work is based on ghost and gluon propagators, which we recently computed on the gauge 
configurations produced by the ETM Collaboration. These are the first to incorporate the effects 
of dynamical charm quarks. In other words, this is an ensemble of sample gauge fields including 
Nf =2+1+1, two light degenerate and two heavy, dynamical quark flavours (see [9, 10] for details 
of the simulations setup) with a twisted-mass fermion action [11, 12]. Combining the propagators 
in Taylor scheme we were able to extract the strong coupling constant, which agrees pretty well 
with the experimental results from T decays [13] and "world average" from PDG [14] at all scales 
from z to Zo. At the same time it became clear that invocation of nonperturbative OPE corrections 
is unequivocally needed to account properly for the lattice data of the coupling in Taylor scheme. 
Within the OPE procedure we expand matrix elements of any non-local operator in terms of local 
operators, organized by their momentum dimensions. Using the sum rules factorization [15, 16], 
these can be expressed as a coefficient to be computed in perturbation (Wilson coefficient) and the 
nonperturbative condensate of a local operator. For such procedure to be consistent it is required 
that (up to a proper renormalization constant for the local operator) the value of the condensate 
should be the same for any Green Function. 

In refs. [17, 18] for Nf=2+\+\ and in refs. [7, 19] for Nf=0 and Nf =2 quark flavours respec- 
tively, we calculate the running coupling by equating numerical results and the OPE prediction. 
The OPE is clearly dominated by the landau gauge gluon condensate A 2 which has been elsewhere 
very much studied (e.g., see [20, 21, 22, 23, 24, 25, 26]). It can be noted that to calculate the 
coupling more reliably we combined the bare lattice propagators in such a manner that the cut-off 
dependence is minimal, which prevents the check for the universality of condensates, as gluon and 
ghost propagators are not to be independently analyzed. However our previous studies in quenched 
approximation suggested that such universality (even including the three-gluon vertex) is indeed 
compatible with the numerical data (Nf=0) [20, 27, 28, 6, 29, 7]. This has been also supported by 
the calculations with Nf-2 dynamical flavours [19, 30]. Now we will address this question in the 
Nf =2+1+1 case, to clarify the nature itself of the unavoidable nonperturbative corrections and the 
involved condensates. 

2. Basics and Simulations 

The Taylor coupling can be obtained directly from the gluon and ghost dressing functions it 
follows [7], 

MP 2 ) = 8 -^ F\p 2 ,a- 1 ) G(p 2 ,a- 1 ) . (2.1) 
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Substituting the dressing functions by their OPE expressions we are left with 

a(p 2 )\™ g 2 (M 2 )(A 2 >R^ 
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where * is known to four loops. 

We refer to the upcoming paper [3 1 ] for the details of this calculation. 
The ghost and gluon propagators are defined as follows: 
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where A" is the gauge field and M = d^D"' 7 is the Fadeev-Popov operator. The lattice setup is 
described in [18] and summarized in Tab. 1. The 0(4)-breaking lattice artefacts have been cured 
by the so-called H (4) -extrapolation procedure [3, 32]. These estimates for bare gluon and ghost 
dressing functions appear plotted in Fig. 1. 
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Table 1: Lattice setup parameters for the ensembles we used in this paper: fc cr j t is the critical value for 
the standard hopping parameter for the bare untwisted mass; /i; stands for the twisted mass for the two 
degenerated light quarks, while ji a and jig define the heavy quarks twisted masses; the last column indicates 
the number of gauge field configurations exploited. This implies that the sttange quark mass is roughly set 
to 95 MeV and the charm one to 1 .5 1 GeV (in MS at 2 GeV), while degenerate light quark masses range 
from 20 to 50 MeV (The lightest pseudoscalar masses approximately range from 270 to 510 MeV). 



3. Testing the OPE 

The gluon and ghost propagator lattice data shown in the upper plots of Fig. 1 are still affected 
by (9(4)-invariant lattice artefacts. As we did in ref. [17, 18, 30] for the Taylor coupling and the 
vector quark propagator, the gluon propagator renormalization constant Z\ att in MOM scheme can 
be written as 



zLatt/ 2 -1 



,a 



2„2 



G(p ,a ) + c a 2pia p 



(3.1) 



which corresponds to the lattice gluon dressing function. In r.h.s., the artefacts-free gluon dressing 
is noted as G. It can be shown that we can rewrite Eq. (3. 1) as: 
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Figure 1: Bare gluon and ghost propagators as they result from lattice computations before (top) and after 
being cured of O (4) -invariant lattice artefacts and properly rescaled to appear superimposed (bottom). The 
blue solid line correspond to the best fit with the OPE formula. 



with 



-NP 



M' 



a(M z ) = G, 



MOM 




r 



,a(p 2 ) 




(A 2 ) R y- 
4 (A? -1);* 



(3.3) 



where Wilson coefficient C w ( ^2,a(jJ, 2 ) J, in the appropriate renormalization scheme, is known 
at the four-loop order with the help of the perturbative expansions for the gluon propagator MOM 
anomalous dimension and Taylor scheme beta function [33]. G NP is expected to describe properly 
the nonperturbative running of the gluon dressing from lower bound of p = 4.5 GeV, below which 
the lattice data deviate from the nonperturbative prediction only including the OPE leading power 
correction, to the upper bound of a 2 p 2 ~ 3.5, above where the impact of higher-order lattice arte- 
facts cannot be neglected any longer. Now we can perform a linear fit and determine the overall 
factor, z(/x,a _1 ), for each data set and the coefficients c a 2 P 2, which should be almost the same as 
they depend only logarithmically on the lattice spacing a(j3) by dimensional arguments (see tab. 2). 
This way we can remove the 0(4)-invariant artefacts away from the lattice data for all available 
momenta in the aforementioned range, by applying Eq. (3. 1), to be left with the continuum nonper- 
turbative running for the gluon dressing; while the overall factors z(n,a~ l ) can be used to rescale 
the data obtained at different /3 such that they all will follow the same curve (see down plots of 
fig. 1). 

As G NP depends on A^r and the gluon condensate, g 2 (A 2 ), their values have to be known 
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Table 2: The overall factors, zQu 2 ,a ), and the coefficients c a 2 P i, for optimal values of g 2 (A 2 ). As dis- 
cussed there, c a i P i is thought not to depend very much on j3 and so is fairly obtained from the results shown 
here. The renormalization point, fi, for the overall factor is fixed at 10 GeV. Statistical errors have been 
computed by applying the Jackknife method. 

prior to the removal of 0(4) -invariant artefacts. In practice, we will take A^g to be known from T 
decays (o^(ml) = 0.334(14) [13], i.e. A^ = 311 MeV for N f = 4) and will search for a value 
of g 2 (A 2 } such that one gets the best matching of the rescaled gluon dressings from the three data 
sets. The details and the intermediate results of this procedure will be explained in [31]. Here we 
just present the results in Tab. 3. 
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g 2 (A 2 ) [GeV 2 ] 


4.7(1.6) 


3.1(1.1) 


4.5(4) 


3.8(1.0) 



Table 3: Results for g 2 (A 2 ) independently obtained from gluon and ghost propagators in this paper and from 
the Taylor coupling analysis in refs. [17, 18]. 

Finally, one can impose the condensates to be the same for both gluon and ghost propagators, 
which leads to g (A ) =4.1 GeV 2 , that appears to be again in a pretty good agreement with the 
previous results in refs. [17, 18]. In the following, we will use the artefacts-free gluon and ghost 
dressing functions, obtained with this last value for the gluon condensate, shown in Fig. 1. 

Now we can directly test the validity of the OPE to describe consistently and accurately dif- 
ferent Green functions by confronting the directly extracted numbers with theoretical predictions: 
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^ 7 W 7 4(JV 2 -1) 

where the l.h.s. is computed after the removal of the artefacts. 

On Eq. (3.4)'s r.h.s., the running with momenta can be easily expressed as a perturbative series 
in terms of the Taylor coupling, being proportional to the gluon condensate, g 2 (A 2 ), which we 
already calculated. As one can see from the Fig. 2, the nonperturbative corrections are totally 
unavoidable for both the gluon (top) and ghost (bottom) cases. Also, the Wilson coefficient derived 
at the order &(a 4 ) in sec. 2 gives a high-quality fit of data (blue solid lines in both plots). Had 
we not included the Wilson coefficient running (i.e., dropped C w from Eq. (3.4)'s r.h.s.), or just 
included its leading logarithm approximation, a flatter slope would be obtained, not describing 
properly the data (blue dotted and green solid lines). A similar situation can be observed in the 
case of Taylor coupling itself, see Fig.3. 
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Figure 2: The left plots shows the difference between the rescaled lattice data for the gluon (top) and ghost 
(bottom) dressing functions and their perturbative estimate, multiplied by the square of momentum, plotted 
versus the perturbative Taylor coupling at the same momentum. For perturbative evaluations, Aj^ is taken 
to be known from T decays. Blue dotted lines correspond to using a leading logarithm approximation for 
Wilson coefficient, and green solid line describes the purely perturbative ansatz. The right plots zoom into 
the region, from 4.5 to 6.3 GeV, where the behavior, predicted by the OPE analysis, is observed. The dotted 
orange line shows a linear fit of data that shows a good agreement with the predicted running. 

Acknowledgments 

We thank for the support the Spanish MICINN FPA201 1-23781 and "Junta de Andalucia" 
P07FQM02962 research projects, and the CC-IN2P3 (CNRS-Lyon), IDRIS (CNRS-Orsay), 
TGCC (Bruyeres-Le-Chatel) and CINES (Montpellier) as well as for GENCI support under Grant 
052271. Speaker is supported by the P2IO LAB EX initiative of France. 

References 



[1] B. Alles, D. Henty, H. Panagopoulos, C. Parrinello, C. Pittori, et al. a s from the nonperturbatively 
renormalised lattice three gluon vertex. Nucl.Phys., B502:325-342, 1997. 

[2] Philippe Boucaud, J. P. Leroy, J. Micheli, O. Pene, and C. Roiesnel. Lattice calculation of alpha(s) in 
momentum scheme. JHEP, 9810:017, 1998. 

[3] D. Becirevic, Ph. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodriguez-Quintero, and C. Roiesnel. 
Asymptotic behaviour of the gluon propagator from lattice QCD. Phys. Rev., D60:094509, 1999. 

[4] D. Becirevic, Ph. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodriguez-Quintero, and C. Roiesnel. 
Asymptotic scaling of the gluon propagator on the lattice. Phys. Rev., D61:l 14508, 2000. 



OPE: ghosts, gluons and CC S 



Konstantin Petrov 



T 



T 



T 



Qh 



0.5 



• p=2.10 
. P=1.95 
P=1.90 

peit+OPE (A MS =320 MeV; g'<A'>=3.85 GeV') 

__ pert (A MS =320 MeV) 




3 4 

p (GeV) 



Figure 3: The Taylor coupling computed with the artefact-free lattice dressing functions and the best-fit 
parameters for A^r and g 2 (A 2 ) (solid blue line). For the sake of comparison, the four loop perturbative 
prediction, evaluated for the same A^j^, is also plotted with a dashed blue line. 

[5] Lorenz von Smekal, Reinhard Alkofer, and Andreas Hauck. The Infrared behavior of gluon and ghost 
propagators in Landau gauge QCD. Phys. Rev. Lett., 79:3591-3594, 1997. 

[6] Philippe Boucaud, J. P. Leroy, A. Le Yaouanc, A.Y. Lokhov, J. Micheli, et al. Non-perturbative power 
corrections to ghost and gluon propagators. JHEP, 0601:037, 2006. 

[7] Philippe Boucaud, F. De Soto, J. P. Leroy, A. Le Yaouanc, J. Micheli, O. Pene, and 

J. Rodriguez-Quintero. Ghost-gluon running coupling, power corrections and the determination of 
Lambda(MS-bar). Phys.Rev., D79:014508, 2009. 

[8] Ph. Boucaud, J. P. Leroy, A. Le Yaouanc, J. Micheli, O. Pene, et al. The Infrared Behaviour of the 
Pure Yang-Mills Green Functions. Few Body Systems, 2012. 

[9] R. Baron, Ph. Boucaud, J. Carbonell, A. Deuzeman, V. Drach, et al. Light hadrons from lattice QCD 
with light (u,d), strange and charm dynamical quarks. JHEP, 1006:111, 2010. 

[10] R. Baron, B. Blossier, P. Boucaud, J. Carbonell, A. Deuzeman, et al. Light hadrons from Nf=2+1+1 
dynamical twisted mass fermions. PoS, LATTICE2010:123, 2010. 

[11] Roberto Frezzotti, Pietro Antonio Grassi, Stefan Sint, and Peter Weisz. Lattice QCD with a chirally 
twisted mass term. JHEP, 08:058, 2001. 

[12] R. Frezzotti and G C. Rossi. Twisted-mass lattice QCD with mass non-degenerate quarks. Nucl. 
Phys. Proc. Suppl, 128:193-202, 2004. 

[13] Siegfried Bethke, Andre H. Hoang, Stefan Kluth, Jochen Schieck, Iain W. Stewart, et al. Workshop on 
Precision Measurements of a s . 201 1. * Temporary entry *. 



[14] K. Nakamura et al. Review of particle physics. J. Phys., G37:075021, 2010. 



OPE: ghosts, gluons and a s Konstantin Petrov 



[15] Mikhail A. Shifman, A.I. Vainshtein, and Valentin I. Zakharov. QCD and Resonance Physics. Sum 
Rules. Nucl.Phys., B 147:385-447, 1979. 

[16] Mikhail A. Shifman, A.I. Vainshtein, and Valentin I. Zakharov. QCD and Resonance Physics: 
Applications. Nucl.Phys., B147:448-518, 1979. 

[17] B. Blossier, Ph. Boucaud, M. Brinet, F. De Soto, X. Du, et al. Ghost-gluon coupling, power 
corrections and A^j s from lattice QCD with a dynamical charm. Phys.Rev., D85:034503, 2012. 

[18] B. Blossier, Ph. Boucaud, M. Brinet, F. De Soto, X. Du, et al. The Strong running coupling at T and 
Z mass scales from lattice QCD. Phys.Rev.Lett., 108:262002, 2012. 

[19] B. Blossier et al. Ghost-gluon coupling, power corrections and A^g from twisted-mass lattice QCD at 
Nf=2. Phys. Rev., D82:034510, 2010. 

[20] Philippe Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene, and J. Rodriguez-Quintero. 
Consistent OPE description of gluon two point and three point Green function? Phys. Lett., 
B493:3 15-324, 2000. 

[21] F. V Gubarev and Valentin I. Zakharov. On the emerging phenomenology of <(A(a)(mu))**2(min)>. 
Phys. Lett., B50 1:28-36, 2001. 

[22] Kei-Ichi Kondo. Vacuum condensate of mass dimension 2 as the origin of mass gap and quark 
confinement. Phys.Lett., B514:335-345, 2001. 

[23] H. Verschelde, K. Knecht, K. Van Acoleyen, and M. Vanderkelen. The Nonperturbative groundstate 
of QCD and the local composite operator A 2 . Phys.Lett., B516:307-313, 2001. 

[24] D. Dudal, H. Verschelde, and S.P Sorella. The Anomalous dimension of the composite operator A 2 in 
the Landau gauge. Phys.Lett, B555:126-131, 2003. 

[25] Enrique Ruiz Arriola, Patrick Oswald Bowman, and Wojciech Broniowski. Landau-gauge 
condensates from the quark propagator on the lattice. Phys.Rev., D70:097505, 2004. 

[26] David Vercauteren and Henri Verschelde. A Two-component picture of the A 2 condensate with 
instantons. Phys.Lett., B697:70-74, 2011. 

[27] Philippe Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene, and J. Rodriguez-Quintero. Testing 
Landau gauge OPE on the lattice with a gluon condensate. Phys.Rev., D63: 114003, 2001. 

[28] F De Soto and J. Rodriguez-Quintero. Notes on the determination of the Landau gauge OPE for the 
asymmetric three gluon vertex. Phys.Rev., D64:l 14003, 2001. 

[29] Philippe Boucaud, F de Soto, J.P Leroy, A. Le Yaouanc, J. Micheli, et al. Artefacts and A 2 power 
corrections: Revisiting the MOM Zpsi (p**2) and Z(V). Phys.Rev, D74:034505, 2006. 

[30] B. Blossier et al. Renormalisation of quark propagators from twisted-mass lattice QCD at Nf=2. 
Phys. Rev, D83:074506, 2011. 

[31] B. Blossier et al. Testing the OPE Wilson Coefficient for gluon condensate from lattice QCD with 
dynamical charm. Preliminary. 

[32] F de Soto and C. Roiesnel. On the reduction of hypercubic lattice artifacts. JHEP, 0709:007, 2007. 

[33] K.G Chetyrkin and A. Retey. Three loop three linear vertices and four loop similar to MOM beta 
functions in massless QCD. 2000. 



